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NEW INTERPRETATIONS OF SHOCK-ASSOCIATED NOISE 


WITH AND WITHOUT SCREECH 

by U. von Glahn 

Lewis Research Center 
Cleveland, Ohio 


SUMMARY 

Anomalous trends in present convergent nozzle (Mach 1) shock-associated 
noise analyses and predictions, with particular t liphasis on the roles of 
screech and jet temperature, are discussed. Experimentally measured values 
of shock-associated noise are used to reassess data trends, including both 
frequency and sound pressure level. The data used includes model-scale 
nozzles, varying in nominal diameter from 5 cm to 13 cm, and full-scale engine 
nozzles up to 48 cm. All data were obtained at static conditions. From this 
reassessment of the measured data, new empirical methods for the prediction of 
shock-associated noise are developed. Separate procedures are presented for 
screech-free and screech-contaminated shock-associated noise. In the present 
approach, shock-associated noise spectra are developed from considerations that 
include the peak .sound pressure level and its frequency, the low frequency 
sound pressure level slope, and the high frequency sound pressure level slope 
or roll-off; the latter Is .shown to vary with directivity angle. 


INTRODUCTION 

In many gas turbine engines, the core flow exhausts at above-critical pres- 
sures, hence at supersonic velocities. With a convergent exhaust nozzle, also 
referred to as a Mach 1 nozzle in the literature, the jet is imperfectly ex- 
panded. The rapid variation in pressure across the nozzle exit plane, according 
to reference 1, Is accomplished by means of steady compression or expansion 
waves that propagate downstream producing the characteristic system of shock 
cells. These shock cells are responsible for two noise sources, (1) screech 
tones and (2) shock-associated noise. 

Screech consists of a discrete spectrum of harmonically related tones and 
is an acoustic feedback phenomenon between the shock cells and the nozzle exit. 
It is also stat. 1 in reference 1 that screech in hot jets is not usually signi- 
ficant. 

When the system of shock cells is sufficiently extensive (pressure ratio 
significantly above choke), shock-associated noise characterized by a broad- 
band, strongly peaked spectrum is produced. This type of noise has been 
studied experimentally in reference 2. From these data, an empirical model 
of shock-associated noise has been developed in the reference. This work does 
not appear to depend on the stagnation temperature of the jet or the jet mix- 
ing velocity. The empirical model proposed was developed by extending a model 
for screech proposed in reference 3, 


In recent years, a significant increase in the shock-associated noise data 
bank has been obtained, particularly at elevated jet temperatures, whereas many 
of the previous studies were made with cold flow. It has also been demonstrated 
(refs. 4 to 7) that the feedback loop associated with screech is easily broken 
by disrupting the smooth trailing edge of the nozzle by the use of tabs at the 
nozzle exit. Care, however, must be exercised to insure that these tabs do not 
adversely affect the jet mixing noise and the aerodynamics of the flow, and that 
their presence does not cause additional shock- associated noise (ref. 7). 

On the basis of the latest studies (refs. 4 to 9), some of the original 
concepts and proposed trends with flow conditions, as discussed later, no longer 
appear to be valid. In the present study, experimentally measured values of 
shock-associated noise are used '•o reassess data trends. Including both fre- 
quency and sound pressure level, The data bank used includes model-scale 
nozzles, varying in nominal diameter from 5 cm to 13 cm, and full-scale engine 
data obtained with nozzles up to 48 cm in diameter. All these data were meas- 
ured at static conditions. From this reassessment and interpretation of the 
measured data, new empirical methods for the prediction of shock-associated 
noise are developed. Separate procedures are presented for screech-free and 
screech-contaminated shock-associated noise. In the present approach, shock- 
associated noise spectra developed from considerations that include the peak 
sound pressure level and its frequency, the low frequency sound pressure level 
slope, and the high frequency sound pressure level slope or roll-off. 


BACKGROUND 

With cold flow, screech tones are easily identified in jet noise spectra. 

A typical example of a broadband screech tone is shown in figure 1 in which 
the sound pressure level is shown as a function of frequency at a 40° 
directivity angle. The data shown are for a 9.5 cm diameter nozzle operated 
with a nominal pressure ratio of 2.75 (ref. 10). Screech tones are generally 
omni-directional so that screech- tone data obtained at one directivity angle 
are similar to those at different angles. In reference 8, narrowband analysis 
(20 Hz) at 9 * 45° were performed to identify the screech tones for a conver- 
gent (Mach 1.0) nozzle. The results of the measured screech frequencies were 
presented in terms of Helmholtz number H « f s D/a 0 . The plot, reproduced from 
reference 8, is given in figure 2. Superimposed on the narrowband data are 
the identifiable screech tones observed in the broadband data of reference 10, 
and unpublished NASA data obtained with 7,6- and 10.2-cm diameter convergent 
nozzles. The broadband screech data fall well within the range of the narrow- 
band measurements. 

As has been discussed in the literature (refs. 4-7), screech tones can be 
suppressed by interrupting the feedback mechanism between the shock cells and 
the nozzle exit through the use of one or more tabs at the nozzle lip. In 
both references 4 and 7 p ; ingle rectangular tab was used to suppress screech 
The rectangular tab had dimensions of 0.125 D wide by 0.063 D long. Compari- 
sons of screech-free and screech-contaminate shock-associated noise are illus- 
trated in figure 3, ta» en from reference 7. In general, the tab eliminated 
the screech tone; however, the broadband jet noise can be somewhat by such a 
device (ref. 7). 
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A frequency parameter for the peak sound pressure level, SPLp, for shock- 
associated noise is commonly given in the literature bys 


where U c is equal to a constant * Uj. The frequency, fp, is defined as that 
frequency associated with the peak SPL in the shock-associated spectrum (SPLp) , 
Hereinafter, fp is called peak frequency. Equation (1) suggests that the peak 
frequency, fp, increases with increasing convective velocity (or jet velocity) 
and decreases with increasing values of jet Mach number (indicated by 8) and 
nozzle size, D. The peak frequency also varies inversely with the directivity 
angle, 9. 

A representative equation of the peak frequency for shock associated noise 
in the literature includes the following: 

f p " c ( so )[( 1 + cos <2> 

The coefficient C varies with the investigator, ranging from 0.763 in refer- 
ence 2 to 0.91 in reference 11. 

Examination of equation (2) indicates that f becomes negative when the 
product of 0.7(U*/a (? )cos 0 exceeds -1.0, as it can .in the rear quadrant 
(9 > 90°). The variation of fp with 0 are given in figure 4 for jet veloc- 
ities of 444 and 922 m/sec, a 8 of 0.83, a nozzle diameter of 9.5 cm and 
C * 0.91. At the lower jet velocity of 444 m/sec (fig. 4(a)), the peak fre- 
quency is positive at all angles because 0.7 (Ui/a 0 )cos 9 does not exceed 
-1,0; the limiting value for the example is -0.91. However, at the higher jet 
velocity of 922 m/sec, this term exceeds -1.0 at a 9 of about 122° (fig. 
4(b)); hence, peak frequencies would be negative at angles greater than 122°. 
This appears unreasonable from physical considerations. 

In order to overcome this deficiency in the formulation of the equation 
for the peak Frequency, an additional term was included in reference 12 where 
fp is given by: 

-1 


( 3 ) 


with C “ 0.7 instead of 0,91 as In equation (2). The variation of f with 
9 is also shown in figure 4 for the same conditions as those used witlvequa- 
tion (2). In general, equations (2) and (3) yield the same trends of f„ with 
9 in the front quadrant (9 < 90°^, however, the computed f p values will be 
somewhat higher using equation ( < than those obtained with equation (2). 
Basically these trends show an increase in f p with increasing 9. In the 
rear quadrant, the peak frequencies computed with equation (2) will be consid- 
erably lower than those obtained using equation (3). With equation (3) a 
reversal in the curve occurs as 9 approaches 180°. This reversal as a func- 
tion of Uj/a 0 . A typical inflection point and curve reversal is shown in 
figure 4(b). Such a characteristic of the peak frequency variation with 
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directivity angle does not appear to be supported by the available experi- 
mental evidence. 

In the preceding formulations for r p , no distinction was made between 
screech-fr»e and screech-contaminated shock associated noise. As will be 
shown later, differences do exist. 


PEAK SHOCK-ASSOCIATED SOUND PRESSURE 
LEVEL FREQUENCY 

As has been shown in the literature (e.g., ref. 7) the presence of a 
screech mode in the spectrum is easily observed with cold flow and low 
3-values but is also present, as will be shown, with heated flow and/or at 
high 3-values. As discussed earlier, a tab at the nozzle lip will eliminate 
the screech and its affect on shock-associated spectra. In the following 
sections, the relationships of the f p with 0 will be developed for both 
screech-free and screech-contaminated shock associated noise. The f p values 
indicated will be for the shock spectra, not the screech tone. 

Screech-free fp. - Values of f p taken from the screech-free model data 
of reference 4 are shown in figure 5 as a function of Uj/ED for several 
directivity angles, 0. The data shown include variations in Uj/a D ranging 
from 0.9 to 1.8. It is seen that f p decreases with decreasing angles for 
constant values of Uj/ED. The line drawn through the data at 0 ■ 90° estab- 
lished a slope with C ■ 0.8. The slopes of the remaining lines shown vary 
from the 0 * 90° line by (1 + cos 0). Good agreement is apparent between the 
measured data and the lines generated by using the following equation for 
determining the f p for screech-free, broadband shock-associated noise: 

f p » 0.8 (Uj/ED) (1 + cos 8) -1 (4) 

A comparison of the f p versus 0 curve resulting from the use of equa- 
tion (4) with those previously shown in figure 4(b) is shown in figure 6. The 
same curve using equation (4) would occur at all Uj/ED values. Thus, the 
shape variations produced by the use of equations (2) and (3) have been elim- 
inated. 

Screech-contaminated f p. - Values of f p as a function of Uj/SD at 
0 - 900 are shown in figure 7(a). The data shown are taken from several ref- 
erences, as indicated in the figure. The solid line is that for screech-free 
f p values taken from figure 5. It is apparent that, for equal values of 
Uj/SD, the screech-contaminated f p are significantly less than those without 
screech. A curve faired through the screech-contaminated data suggests a 
constant of C « 0.6 compared to the value of C * 0.8 for the screech-free 
data (eq. (4)). At directivity angles other than 90°, the trends between the 
screech-free and screech-contaminated f p values are similar. For example, a 
comparison of f p values at 0 * 60° is shown in figure 7(b). The dashed 
line through the data is obtained using the constant of C « 0.6 determined 
from the data at 6 » 90° and using the relationship (1 + cos 0)” 1 as in 
equation (4). The data in figures 7(a) and (b) are for model-scale nozzles. 

In figure 7(c), similar data are shown for full-scale engine nozzles (ref. 12) 
at a nominal 0 of 50°. On the basis of the data shown in figure 7, the f p 
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for screech-contaminated, shock- assoc laced noise can be estimated with the 
following equation: 

f p - O.6(bVj/0D)<l + cos 9)" 1 (5) 

PEAK SHOCK-ASSOCIATED SOUND PRESSURE LEVEL 
Screech-Free SP! V 

r 

In the development of the present analysis. It is assumed that the varia- 
tions of both the peak sound pressure levels, SPLp, and the overall sound 
pressure levels, OASPL, with the jet Mach number parameter, 8, are similar. 

The variation of OASPL with 8 given in reference 9 for the case of cold flow 
and a directivity angle, 9, of 45° is shown in figure 8. From these measured 
data, an empirical relationship was obtained yielding the curve shown in the 
figure. This empirical relationship is expressed by: 


oaspl q „ 45 o 


- 118.5 + 20 log 8-10 log 



( 6 ) 


Also indicated in the Figure by the two dash lines are the slopes (not 
absolute values) of simple approximations for the variation of OASPL (hence- 
SPLp) with 8. In the region 0.4 < 8 < 0.9 the variation of OASPL with $ 
is approximated by 8 . In the region of 6 > 0.9, the approximation is given 
by 82, These exponents of 8 have been used by some investigators to 
describe the variation of SPL p with 8- It is further assumed that with 
heated jet flows, the trend of the OASPL with 8 is also given by equation (6), 
although the constant, 118.5, in equation (6) will change with jet temperature. 


In terms of SPLp, an effective sound pressure level parameter is used 
herein designated SPL*. This latter term is defined as: 

r 

SPLj ~ SPL p - 20 log D/R - F^S) (7) 

where F^(B) is given by: 

F^S) * 20 

The variation of SPL p with 9 is shown in figure 9, the data being 
taken from reference 4. Tne curves shown in figure 9 are based on a simple 
10 log (1 + cos 6) relationship and fit the data well. This directivity term 
is now incorporated into equation (7) yielding: 

, + 0.105 
< 8> 3 * 6 


SPLj - SPLp - 20 log D/R- 20 log 8 + 10 log 


j-lQ log(l + cos 0) 



( 9 ) 


mmmmm 
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TIiub far, the analysis has dealt with cold Jet flows* When the measured 
SPLp values of shock-associated noise are plotted in terms of equation (9) as 
a function of jet static-to-ambient temperature ratio, Tj/T 0 , the SPLp values 
first increase with increasing temperature ratio and then decrease (fig. 10). 
Also the data show a dependency on the jet Mach number parameter, 8, The data 
can be fitted empirically by the curve shown in figure 11 using the following 
relationship: 



153.3 - 15 log F 2 (p) - 10 log 


1 + 


0.7 


(F 2 (8)) 


1.5 


(10) 


where 

F 2 <8) * (Tj/T 0 )/8(1 + O.3750 3 ) (11) 


by 


By combining equations (9) to (11), the value of SPLp can ba calculated 


SPLp « 153.3 - 15 log 


V*. 


8(1 + 0.375B 3 ) 


- 10 log 


1 + 


r 6(l + 


°- 7 ( V* 


„ 4.5 
0.3758 3 ) \ 

/T 0 ) 


+ 20 log D/R + 20 log 8 - 10 log 1 1 + 

<&)' 


Fi + -JELtiaJ 

L (8) 3 - 67 J 


+ 10 log(l + COS 0) 


( 12 ) 


Screech-Contaminated SPLp 

The analysis procedure for screech-contaminated shock-associated SPLp 
is identical to that used for screech-free shock noise in the preceding 
section. In figure 12 is shown the variation of SPLp with the previously 
introduced parameter Tj/T 0 /P(l + 0.375P 3 ). The data shown include both 
model-scale results and those from two engine tests (Orenda engine used in 
the F-86 aircraft (ref. 12) and the J-85 engine used on the Aerotrain 
(ref. 12). The data show the same trend with temperature ratio as that with 
the screech-free data in the preceding section. The solid curve shown in the 
figure is given by: 


SPLj = 154.3 - 15 log F 2 (3) - 10 log 1 + ^7-? (13) 

L (F 2 (8)) 4 * 5 J 

where F 2 (g) is given by equation (11). This equation is similar to equa- 
tion (10) except for the constants. Also shown for comparison in figure 12 
by the dotted line is the curve for screech-free shock-associated SPL?. A 
comparison of the two curves shown indicates that the screech-contaminated 


shock noise is less than that for screech-free shock noise at temperature 
parameter values less than about 1.55. At temperature parameter values 
greater than 1.55 the opposite trend occurs. 


On the basis of the relationships shown in figure 12, the following equa- 
tion can be used to estimate screech-contaminated shock-associated peak sound 
pressure levels $ 


SPLp « 154.3 - 15 log 


V T o 


3(1 + 0.3758 3 ) 


10 log 


(■ 


1 + 2.7 (£O±-0 i 375Bfi 


Tj/T 0 


\4.5 


+ 20 log D/R + 20 log g- 10 log 1+ ~rrr + 10 log(l + cos 0) 

L <e) 3 - 67 J 

(14) 

SHOCK-ASSOCIATED NOISE SPECTRAL SHAPE CHARACTERISTICS 

On the basis of the data referenced herein, the spectral shape character- 
istics are substantially the same for both screech-free and screech-contaminated 
spectra. Consequently, the same empirical relationships are assumed herein to 
apply to both spectra. 

A representative screech-free shock-associated noise spectrum (ref. 4) is 
shown in figure 13. The spectrum shown has the jet mixing noise subtracted 
from the total measured SPL values. The following sections describe the 
spectral characteristics in two regions: (1) that portion of the spectrum in 

which f < f p , and (2) that portion of the spectrum in which f > fp. 

Lower frequency forward spectrum, f < f p . - In general, the data indicate 
that the shock-associate3 SPL in the lower frequency forward portion of the 
spectrum increases with frequency at a rate of about 21 dB/octave band. This 
rate appears to be independent of all the variables included in the data bank. 

Higher frequency aft spectrum, f > fp . - In the higher frequency aft por- 
tion of the spectrum, the shock-associated SPL decreases with increasing 
values of frequency. The variation of SPL with frequency, however, appears 
to be a function of the directivity angle. A typical variation of the change 
in the SPL/f slope with 0 is shown in figure 14. It is apparent that the 
slope of SPL/f has decreased with a change in 0 from 45° to 90°. Although 
the data shown in figure 14 are for screech-free shock-associated noise, simi- 
lar results were obtained with screech-contaminated shock noise. This is to 
be expected because introduction of the tab that eliminated screech in these 
data did not affect the higher frequency portion of shock-associated noise 
spectrum (ref. 7). 

Examination of the higher frequency aft portion of the shock-associated 
spectral data indicated a temperature dependency of the aft SPL slope. 

Typical aft SPL data (ref. 10) in terms of dB/octave band slope as a function 
of temperature ratio, Tj/T Q for three representative directivity angle (0 * 20°, 
60°, and 90°) are shown in figure 15. It is apparent from these data that the 
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aft spectral slope becomes steeper, i.e., greater negative values of dB/octave 
band with increasing temper ature ratio. The dB/octave band varies with 
(Tj/T o )0*25 as shown by the solid curves in the figure. The aft spectral 
slope data from reference 10 are replotted in figure 16 in terms of (dB/octave 
band)/(Tj/T o )0*25 as a function of Tj/T 0 . On the basis of this data plot, 
an equation was developed to correlate this term with directivity angle. The 
correlation equation is given by: 

(■IB/octav. band) . _ 2 n l COB £\ _ 0-095(1 + cob) 5 (15) 

(Tj/To) 0 ' 25 ' 2 > 

The solid curves in figure 16 were obtained using equation (15). In general, 
the comparison between the calculated curves and the data are good. This is 
further shown in figure 17 where the data are plotted in terms of (dB/octave 
band) /(Tj/T D )®*25 a s a function of directivity angle, 0. Shown for the data 
at each angle is a mean value and the range of the data shown in figure 16. 
Also shown is the curve calculated using equation (15). Although the data 
deviate somewhat from the calculated curve at 0 * 40° and 60°, this may be 
the result of attempting to establish a precise data slope from the experi- 
mental data. 


COMPARISONS OF MEASURED AND CALCULATED SPECTRA 

In the following sections, representative measured spectra at selected 
directivity angles from 40° to 90° will be compared with those synthesized by 
the techniques discussed herein. The synthesized spectra are developed from 
the following procedure: 

(1) The SPLp and fp are calculated from equations (12) and (4), re- 
spectively for screech-free shock associated noise or from equations (14) and 
(5), respectively for screech contaminated shock-associated noise. These cal- 
culations yield the apex of the spectral shape. 

(2) The low frequency portion of the spectral shape is obtained using a 
constant 21 dB/octave band slope through the previously calculated spectral 
apex (SPLp and fp). 

(3) The aft spectral roll-off is calculated using equation (15). The 
roll-off is initiated at the spectral apex (SPL p and f p ). 

Based on this procedure, both screech-free and screech-contaminated, 
synthesized shock-associated spectra are compared with data from references 4 
and 10 in figures 18 to 20. The data shown have the jet mixing noise subtrac- 
ted from the total measured SPL values. 

Screech-free spectra . - In figure 18 is shown a comparison of the measured 
and calculated spectra for selected flow conditions. The data are taken from 
reference 4. The spectra shown are for a g * 0.94 and jet velocities of 341, 
538, and 601 m/s. The corresponding Tj/T 0 ratios are 0.73, 1.82, and 2.27, 
respectively. The calculated spectra are for 0 of 45°, 60°, and 90°, with 
data shown where available. It is apparent that, in general, good agreement 
exists between the measured and calculated spectra. 
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In figure 19, a spectral comparison is shown for a B ■ 1.34, a Jet vel- 
ocity of 587 n/s and a Tj/T 0 ratio of 1.82. Good agreement between the 
measured (ref. 4) and calculated spectra is again evident. 

Screech-contaminated spectra . - In figure 20 is shown a comparison of the 
measured and calculated screech-contaminated spectra for selected flow condi- 
tions. The spectra shown are for 8 * 0.83 and jet velocities of 337, 540, 
and 746 m/s. The Tj/T 0 ratios are 0.72, 1.52, and 2.96, respectively. The 
data shown are from reference 10. The agreement between the measured and cal- 
culated values is again as good as that for the screech-free data. 


CONCLUDING REMARKS 

On the basis of the present analysis of available shock-associated noise, 
it appears that, unless screech tabs were used, the measured data were contam- 
inated by screech. The data also Indicate that nozzle size (scale) is not a 
factor as to whether or not screech contaminates the shock-associated noise 
since the noise measurements included nozzles ranging from 5 to 48 cm in 
diameter. The data analysis also appears to indicate that screech persists 
even at elevated jet static temperatures and high pressure ratios. At these 
latter conditions screech is almost unnoticeable. In all cases where screech 
appears to influence the shock-associated noise, the spectra peak at a lower 
frequency, by 1/3- to 2/3-octave band, than those that are screech-free. Fin- 
ally, it is apparent from this study that the rigs used in the several experi- 
mental programs were not a factor in the shock-associated data. 

It should be pointed out that, to the best of the author's knowledge, all 
the nozzles included herein have smooth contours at the exit plane. It is not 
known, at this point, if a full-scale variable area nozzle consisting of over- 
lapping metal leaves or "feathers" would constitute the equivalent of tabs on 
a screech-free nozzle. The Irregularities caused by the leaves or "feathers" 
on the nozzle exit could help to eliminate the feedback loop between the shock 
c^lls and the nozzle lip. The experiments performed with various screech- 
elimination devices in reference 6 would tend to support the assumption that 
such nozzles are screech-free. 


SUMMARY OF RESULTS 

Empirical correlations have been developed that can be used to synthesize 
and predict both screech-free and screech-contaminated shock-associated spectra. 
The predicted spectra are in good agreement with the measured data and reflect 
trends implicit in the data, including the following: 

1. The peak shock-associated sound pressure l evel SPLp , varies with jet 
Mach number parameter 8, where 8 is defined by - i. For jet Mach numbers, 

Mj, less than about 1.4, the SPL p variation is approximated by 8^ whereas 
for Mj values greater than about 1.4, the variation is approximated by 8^. 

An equation was developed herein to provide for the continuous variation of 
SPLp with 8. 
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2. The peak sound pressure level, SPL p , varies with jet static tempera** 
ture. At low jet static temperatures, the SPL p increases with the tempera** 
ture; however, at high jet static temperatures the SPLp decreases with 
temperature. An equation was developed to reflect this trend. 

3. The frequency associated with SPL p is greater by a factor of 1.33 
for screech**free cases than those for screech-contaminated cases. 


4. The spectral characteristics at frequencies less than that at the 
SPLp value vary with an apparent dB/octave band flops of 21. At frequencies 
greater than the frequency at SPL p , the spectral characteristics are a func- 
tion of jet static temperature and directivity angle, 6. An equation was 
developed to reflect this trend in the data. 

5. The variation of both the SPLp and its associated frequency with 
directivity angle, 0, is given by (1 + cos 0). 


SYMBOL LIST 


(All symbols are in SI units except as noted.) 
a 0 ambient speed of vo ;nd 

C constant 

D nozzle diameter 

F,Pj,F 2 function 
f 1/3-octave band frequency 

fp peak sound pressure level frequency 

f g screech frequency 

H Helmholtz number, f g D/a 0 

Mj jet Mach number 

OASPL overall sound pressure level, dB re 20pN/ro^ 

R distance from nozzle exit to microphone locations 

SPL sound pressure level, dB re 20pN/m2 

SPL peak sound pressure level, dB re 20pN/m^ 

r 

SPL p sound pressure level parameter 

Tj jet static temperature 

T 0 ambient temperature 
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convective velocity, con.tent * Uj 


Jet velocity 


Mach number parameter. 



directivity angle measured from inlet 
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Figure Z - Screech tones for convergent (Mach 1) 
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Figure 8. - Overall sound pressure level variation with 0. 
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